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Abstract 

A canonical formalism is presented which allows for investigations of quantum radi- 
ation induced by localized, smooth disturbances of classical background fields by means 
of a perturbation theory approach. For massless, non-selfinteracting quantum fields at 
zero temperature we demonstrate that the low-energy part of the spectrum of created 
particles exhibits a non-thermal character. Applied to QED in varying dielectrics the 
response theory approach facilitates to study two distinct processes contributing to the 
production of photons: the squeezing effect due to space-time varying properties of the 
medium and of the velocity effect due to its motion. The generalization of this approach 
to finite temperatures as well as the relation to sonoluminescence is indicated. 
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1 Introduction 

Perhaps one of the most fascinating aspects of the static Casimir effect |[[ p] is that it demon- 
strates the highly non-trivial nature of the vacuum state in quantum field theory (see Refs. 
||. [J, Hi for reviews). During the last decade the dynamical Casimir effect, i.e. the production 
of particles due to the presence of space-time varying external conditions has become subject 
of an increasing number of studies especially devoted to particle creation in gravitational back- 
grounds (see for a review |5|, §), due to dynamical boundary conditions (see Refs. in @) and 
photon production in moving dielectrics (see e.g. P, |9], |10[) to mention only a very few of 
prominent fields of research. 

The phenomenon of quantum radiation (QR) can be interpreted as the response of the vac- 
uum of quantized fields due to their interactions with external conditions. One may think of 
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the latter as some space-time dependent background fields (classical field configuration) or as 
time-dependent, geometrical boundary conditions superimposed on the quantized fields under 
consideration (e.g. moving mirrors or dielectrics). 

In a recent investigation 0] we propose a canonical approach to the dynamical Casimir effect 
in the presence of boundary conditions which facilitates to calculate the number of generated 
particles by means of response theory. In this article we are going to apply this approach to the 
situation of quantum fields in arbitrary space-time dependent background fields lying particular 
emphasis on the situation of QED in moving and space-time varying dielectrics. Under rather 
general conditions the low-energy behaviour of the spectrum of produced particles can be de- 
duced. We also indicate that the formalism is most suitable for incorporating finite-temperature 
effects as well. This paper is organized as follows: In section 2 we deduce some features of the 
energy spectrum of the dynamical Casimir effect in smooth background fields based on the re- 
sponse theory approach. In section 3 the formalism is applied to QED in dynamical dielectrics 
lying particular emphasis on the problem of field quantization. An interaction Hamiltonian is 
derived which describes particle production due to two distinct generic prossess: 1. squeezing 
of the QED ground state due to localized, space-time varying dielectric properties inside the 
medium; 2. velocity effect due to motion of the medium. Section 4 examines the squeezing ef- 
fect where a general expression for the total radiated energy will be derived. Section 5 discusses 
the velocity effect qualitatively to underline the importance of realistic ansatze for the velocity 
field and how they may affect the energy-spectrum of the produced particles. We will close 
with some discussions devoted to the following issues: generalization of the formalism to finite 
temperatures, application to media with nontrivial permeability function and an indication of 
how to estimate the contribution of QR to the phenomenon of sonoluminescence. 



2 Energy spectrum of quantum radiation 

In this section we present the general formalism for the calculation of the number of particles 
produced by the dynamical Casimir effect. We consider a set of quantized fields 0,4 interacting 
with a classical background described by some smooth, space-time dependent functions AKf^. 
The rather general system under consideration may be characterized by the four suppositions: 

I : The 9^ are supposed to be massless and non-selfinteracting bosonic fields . 

II : The field configuration is considered at zero temperature 

and initially in the vacuum state. 

III : A proper definition of particles is provided by the unperturbed Hamiltonian Hq . 

IV : The interaction Hamiltonian Hi is determined by small, space-time localized and 

smooth functions AK^ . 
These four assumptions allow to deduce the low-energy behaviour of the spectrum e(u) of the 



2 



produced particles. We shall prove that for small freqencies u the energy spectrum e(u) will 
be proportional to u 4 in contrast to the thermal spectrum which behaves as u 3 . 



At first we specify a class of Lagrangians compatible with the four assumptions stated above. 
Due to supposition I we assume £ as a bilinearform of first order derivatives of the fields 0^ 
interacting with background fields AK AB : 



C 



\ d»&A (K* B + AK AB (r, t)) d»Q B (1) 



with AK AB {f,t) -> for | r | — > oo and for t — > ±oo (according to IV). In Eq. flU) and in the 
following formulae we make use of the summation convention and declare that one has to sum 
over all indices (A, /i, i etc.) that do not occour at both sides of the equation. 
This leads to the total Hamiltonian : 

H(t) = H + H x (t) = J d 3 r (Ho + Hi(t)) (2) 

with 7ii(r, t) — > for \ f | — > oo and for t — > ±oo (according to IV) together with the undisturbed 
Hamiltonian density 



1^ 1 
2 ' " 2 l] 



n = -u A T AB n B + -v 1 f B d i e A cPQ B , (3) 



assuming the constants K AB satisfy (possibly after a suitable transformation) K AB = K^ B = 
or other appropriate conditions leading to the constant matrices T AB and V^. 
A proper definition of particles with respect to the undisturbed Hamiltonian Hq (according to 
III) is provided by insertion of the expansions 




ikr 



+ h.c. 



(4) 



for the fields and 



U A {r,t) = V 



-1/2 y^B 



h.cO 



(5) 



for the canonical field momenta. (V denotes the quantization volume.) 
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Supposition III demands the existence of appropriate matrices F- and Gg which allow to 
diagonalize H (i.e. terms like o,^ A ap B and o,tAt/ B generated by the insertion of and 
into (pf) have to cancel) leading to the form 



(6) 



kA 



kA 



The i^o-frequencies LOt A possess a quasi-linear dispersion relation (after f) : 



, M v /.v/.- ; ii;V o{k k) . (?) 

The following calculations are most suitably performed in interaction-representation: 



dA ... j. &4 



(8) 



(9) 



Suppositions II and III imply the following definition of the vacuum - the initial state - as well 
as that of proper particle creation/annihilation operators o^ A /a^ A : 



\ip(t -> -oo)> = |0) with V kA a u |0) = . 
Eq. (P) can be formally integrated with the time-ordering operator T 



(10) 



m -> oo)) 



T 
T 

cx 

E 

n=0 
oo 

E 

n=0 



exp ^— i J dt H\{t) 
exp I —% j d 4 x Hi (x) 



10} 



oo f_* n 



ni 
-i) 



nl 



10) 

dt n ---J dhT [HxiQ ■ ■ ■ Htitt)] |0> 
d A x n --- j d A x\T TCi(x n ) ■ • ■Hi(x 1 ) |0) 



with x = (r,t). The expectation value of the number opertator is given by 



E 

m=l 



ml 



d^xi 



n=l 



nl 



4„/ 



da; 



x_ m ) 7^ \ {x_\ j 



n %a t [Md n )---MdS\ |o> 

4 



(12) 



The vacuum expectation values have to be calculated in the ifo-dynamics and can also be 
deduced from the 2(n + m + l)-point correlation functions (0| Q A i(xi) 4 4 4 ®A v {x p ) |0). Due to 
assumption I the pertubation Hamiltonian density can be expressed as: 



H 1 = Q^(x)d»Q A d»Q B . (13) 

In the interaction representation the time evolution of the fields is now governed by Hq an d 
assumes the form: 



e A (x) = v-wf* f ^^tHb e + h - c -) > ( 14 ) 

k 

with k A = (k,uj^ A ). The interaction term becomes 



ftxHiid = ± . KcKd x 

J rf 4 xg^(x)Fi ? Fi , (a| iCe ^ + h.c.)(a+ D e^ + h.c.) . (15) 

In order to obtain a non-zero matrix-element in (^) at the left and at the right of N^ A at least 
one summation index kg has to be equal to k, which generates an overall factor of order k . 
Assuming the functions Q^(x) to be localized and smooth enough (IV) , all mode summations 
are convergent and the d 4 x- integrals are finite for the limit k — > 0. Under these conditions, we 
can conclude that the number of particles per mode possess the low-momentum behaviour 

{N U ) = 0{k) . (16) 

This enables us to calculate the total radiated energy associated with the total number of 
created particles: 

OO CO 

E = f W£a (N U ) - £ / d 3 k W£A (N U ) =Y.fdu A e A {u A ) =Jdu e{u) . (17) 



A A 



kA 







In the limiting case of small frequencies u simple power counting leads to a non-thermal spectral 
energy density: 

e(w) ~ + 0{u 5 ) . (18) 

This reveals a general feature of the low-energy part of the spectrum of particles created by 
the dynamical Casimir effect due to the interaction with external background fields, assuming 
that the suppositions I - IV , which describe a very general case, are fulfilled. 
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3 QED in dielectrics 



3.1 Equations of motion 

Throughout this article natural units with h = c = eo = po = l will be used. We start with 
the source free Maxwell equations ( i.e. for p = and j — ) : 



VB = , V x E = -B , (19) 

VI) = , V x H = D . (20) 

The first pair ( |i~9l) leads to the definition of the potentials : 

B = -VxA , E=A+V<5> . (21) 

Now we suppose a linear and nondispersive medium in arbitrary motion described by the scalar 
dielectricity function e(r, t) and a velocity field /3(f, t). For the special case of a medium globally 
at rest, i.e. — (1,(3 — 0), where the relations D = eE and H = B hold, the Lagrangian 
leading to the second pair of the Maxwell equations reads: 

C = \{eE 2 -B 2 ) . (22) 

There are different options for the construction of the Lagrangian for the general case of arbitary 
u^(f,t) and e(f,t). The most direct way is to make an ansatz as a linear combination of all 
possible Lorentz-covariant scalars constructed from the tensors at hand: the field strength 
tensor F^, the Lorentz metric g^ v , the four velocity u^, possibly the Levi-Civita pseudo-tensor 
£nupa and the dielectricity function e (Lorentz-scalar). For reasons of simplicity one may neglect 
additional acceleration effects that could be carried by u^ v or e ^ etc. The only form that 
reduces to the special case ( p2|) of a medium at rest and the vacuum Lagrangian £ = ^F^ u F Ufl 
for e — 1 , is provided by 



= ±(E 2 - B 2 ) - - (0Ef - (E - x Bf) . (23) 

This Lagrangian established by means of covariance arguments may be expanded in powers in 

— * — # 

(3 (assuming \(3\ ^ 1) resulting in the non-covariant form 

C = ^(sE 2 -B 2 ) + (s~l)P(ExB) + O0 2 ) . (24) 
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The Euler-Lagrange equations for the scalar potential $ lead to: 

vJj = Vi3 = (25) 

while the equations of motion for A read: 

d dC ■ d dC ,„ t t \ / \ 

o - = A = -t^itt- = V x HI . 26 

So the Euler-Lagrange equations are identical with the second pair of the Maxwell equations, 
if we introduce the electric and magnetic displacement fields D and H according to 

BC -> 

D = ^ = eE + (e-l)Bxp + 0(p 2 ) (27) 
dE 

and 

H = -^ = B + (e-l)Ex(3 + 0((3 2 ) . (28) 
Applying Euler's theorem the Lagrangian ([23|) transforms into the alternative form 

C= l -(ED-BH) , (29) 
which holds to all orders in j3 . 

3.2 Gauge problems 

If we try to carry out the canonical quantization based on the Hamiltonian formalism applied 
to the Lagrangian obtained above, we are faced with some problems arising from the fact, that 
QED is a gauge field theory. (Attempts for quantizing QED in dielectric media date back for 
more than half of a century [JTTJ, |T2[. To solve the quantization problem for general media has 
attracted recent interest in quantum optics, see |13|-||20||.) Due to the fact, that £ contains no 



kinetic term for the scalar field component $ there exists the primary constraint: 

n = ^ = o . (30) 

A secondary constraint arises from the absence of any potential term for the longitudinal 
component of the vector potential An: 

VI1 = V— = V^S = VD = . (31) 

dl dE 

So far these problems also occur in vacuum QED. However, in that case it is possible to remove 
both kinds of field modes simultaneously by choosing the Coulomb gauge VA = which has 
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$ = as immediate consequence. Equivalently, one can choose the temporal gauge $ = 
which directly leads to V 'A = 0. For a detailed discussion of the problem of the quantization 
with constraints we refer e.g. to [2T|. So in vacuum QED it is possible to remove the scalar 



and the longitudinal photons simultaneously. This fails however for QED in media because of 
the non-trivial relation between E and D. As the next step we observe that the Hamiltonian 
derived according to the canonical prescription 

tt(n,I,$) = uA-c = -(ed + bh) -nv$ (32) 

contains the term II V$. Even though VII = 0, dropping this term or setting $ = would 
change the equations of motion, e.g. such that they are no more consistent with the Coulomb 
gauge V 'A = 0. One way to solve the gauge problem is provided by carrying out the procedure 
of Dirac quantization as described e.g. in One could also try to generalize the path integral 
quantization recently applied to dielectrics at rest [20|. Here we proceed in a different way: 



Based on a decomposition of the electromagnetic fields into an appropriatly choosen basis {f a } 
we introduce a set of independent variables {q a } leading to a Lagrangian L(q, q, t) without any 
constraints and perform the usual canonical quantization. It will be necessary to prove that 
the primary and secondary constraints required for the fields can be satisfied by the adequate 
choice for {/ a }- 

For this purpose we fix the gauge by choosing the Coulomb condition: 

VA = . (33) 

As mentioned above, in dielectrics this does not automatically imply $ = 0. Furthermore, from 
the equation of motion V-D = it follows: 

V(eV$) = -V (e A +(e -l)^x(VxI))+ O0 2 ) , (34) 

which possesses non-trivial solutions $ 7^ 0. Requiring the boundary conditions $ — > for 
|r I — > 00 Eq. (|34l) is solvable with the aid of the Green function G(f,r',t) corresponding to 
the Laplace-Beltrami-type operator (VeV) : 



$(f,t) = - / dV G(f,f',t) x 

V (e(r',t) A (r',t) + (e(f',t) - l)(3(r',t) x (V x A(f' ,t))) + O0 2 ) . (35) 



This equation explicitly reveals the nature of the scalar potential $ as a dependent field vari- 
able. Since Eq. (31) represents an equation of motion the action A becomes extremal at 
the hypersurface described by (|35]) . Therefore the correct equations of motion could also 
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be obtained by varying A only on this hypersurface. Accordingly, it is allowed to elimi- 
nate $ in favour of the vector potential A by inserting Eq. (^) into L[A, A, Q,t] arriving 



at L'[A,A,t] = L 



A,A,$[A, A,t],t 



The next important step is to expand the remaining fields into an appropriate basis set exploit- 
ing the gauge condition. The irreducible degrees of freedom of A with VA = were described 
by an expansion into a complete set of real, orthonormal and transverse functions / (r) (see 
appendix): 



A(f,t) = q a (t)f a (f) , q a (t) = J d 6 rA(r,t)f a (r) . (36) 

Insertion of this expansion into ( j35|) allows to express $ completely in terms of the variables q 
and q : 



$(r,t) = q a (t)ip a (f,t) + q a (t) X a(r,t) , (37) 
<Pa(f,t) = -J d 3 r' G(r,r',t)V ((e(f',t) -1)0 (f',t) x (V x ^(f'))) + 0(/3 2 ) ,(38) 

Xo (f, t) = - f d\' G(f, r ', t)V (e(r £)/>"")) + 0(/3 2 ) . (39) 



The scalar potential does no longer appear as a degree of freedom. We can express the physical 
fields: 

E(r,t) = q a {t)f a {r) + (la{t)V Va {r,t) + q a {t)VXa{r,t) , (40) 

B(f,t) = -q a (t)Vxf a (f) . (41) 

Insertion of these decompositions into the Lagrangian Eq. (E3) leads to 



L(q, q,t) = Jd 3 rC = J d 3 r (Co +A(t) 



= 7}( £ °°<il - ^l<ll) + q a X ab (t)q b + q a y a b(t)qb + q a Z ab (t)q b . (42) 

Now we have a Lagrangian expressed by an independent set of variables q, that contains no 
constraints. The Euler-Lagrange equations 
d dL dL 



dt dq a dq a 



(43) 



lead to the Maxwell equation D = V x H after multiplication with f a (r ) and reidentification 
of the fields. The other equations - especially V-D = - are automatically satisfied. Now we 
can derive the canonical momenta p and their field representations II (see appendix): 

n=p a f a = —f a = (eE + (e-l)BxP + O0 2 )) ± = D . (44) 
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It turns out that the canonical momentum is equal to D also in this general case. This is not 
obvious because the terms including <3> contain also q ; but performing the calculations we see 
that the contribution from q,t) to p vanishes due to VD = (see appendix): 
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6 A, 



-L 



A,A,$[A,A,],t 



f a —L(q,q,t) = -^L \E,B,t 
dq a 5E 1 



(45) 



The Hamiltonian now becomes: 



H(p, q, t) = p a q a - L — J d 3 r H(p, q, t) 



(46) 



with 



H{p,q,t) 



respectively 



1 ( 1 



2 \e(f,t) 



n 2 + (V x A) 2 + 



e(r,t)-l 
e{f,t) 



0(f, t) (n x (V x A)) + O0 2 ) (47) 



H(p,q,t) = ~(ED + M) 



(48 



We have thus derived the Hamiltonian which indeed has the same form as the one obtained 
in the naive way (0), except for the term II V$. Now $ and Aii are no longer degrees of 
freedom and the gauge condition is automatically satisfied. As a consequence the term II V$ 
disappears. 

3.3 Quantization 

Based on the Lagrangian L(q,q,t) obtained above it is possible to carry out the procedure of 
the canonical quantization via defining the following commutation relations: 



[Qa(t),qb(t)} = \pa(t),p b (t)} = , [q a (t),Pb(t)} =i5 ab . 



(49) 



Multiplication with the eigenfunctions {f a } leads to the commutation relations for the fields: 



[i*(r,t),^(f',t)] = [n l (f,t),fp (f',t)} = 



(50) 



and 



[A\r } t) } il j (r\t)} =i5 i l(f-f') 



(51) 
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Eq. (|5T| ) can be used to obtain a medium-independent (i.e. for arbitary s(r,t) and (3(f,t)) 
commutation relation between the fields B = — V x A and D = fl (see e.g. : 

[^(r^Z^f',*)] =ie <J '*Sl fe <y(F-F / ) . (52) 

In the quantization prescription based on the decomposition presented above the Coulomb 
gauge condition (and also V-D = 0) appears as an operator identity. This leads to 

V |V>),n : m(Vl) n m = , (53) 

i.e., the fixing of the gauge removes all longitudinal photons. But for the scalar photons this is 
not the case for any variation of e(f,t) or /3(f,t) : 

(0|l> 2 (f,t)|0>^0 . (54) 



3.4 Interaction Hamiltonian 

Now we are in the position to apply the formalism presented in section 2. After separating the 
undisturbed Hq , the interaction Hamiltonian can be cast into the form: 

*• = \(wa - i) *' + eJ wjr i?ir ' t] (* x (v x X) ) + • (55) 

There are two terms giving rise to disturbances of the vacuum and thus contributing to the 
production of photons: the first term due to the change of the dielectric function e(f,t) - the 
squeezing term - and the second one related to /3(r, t) - the velocity term. 
The assumption of localized perturbations as stated in section 2 requires e(r,t) and 
P(r, t) — > for |r | — > oo and for t — > ±oo. 



4 Squeezing effect 

As indicated in section 2, the general form for the particle production due to the disturbance 
Hi turns out to be: 



= Jd 4 xJ d 4 x' (0| H^N^Hiix) |0> + O(Hf) = + 0(H*) . (56) 

Now we shall evaluate this expression for (3 = 0, neglecting the effects of any matter flow, i.e., 
for the pure squeezing contribution: 

Mr, t) = \ (-^ - j^j n (r, t) = t) n(f, t) . (57) 
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The mode expansion in the ifo-dynamics leads to: 



n 



fcA 



with k = (k, oj^) = (k,k/ y/e^Z) and ke^ x = . Inserting this expansion into H 1 we arrive at : 



6- Oj- 6^ 6 l - * - — 2)— 
Al /ci Ai k2 X'2 &2 A2 



felAl fc2A2 



x^mf f "^"iJ^L/^/^- W -(59) 

£3-^3 fc4A4 

Performing all /cA-summations we would obtain a distribution-like response function that is 
related to the 6-point correlation function. Here we proceed in a different way and carry out 
the rf 4 x-integration which yields a more compact form of the equation in terms of Fourier 
transformation JF: 

T : f(x) - /(£) = / d 4 x f(x) e** . (60) 
Accordingly, we obtain 



^ = f f f f — rfe + fc 2 )efe + *u) x 

fclAl /C2A2 &3A3 fc4A4 

(% 1 A 1 %A 2 )(%3A3%4A 4 ) (°l«fc 1 A 1 «fc 2 A2^A«£ A 3^ 4A4 1°) ( 61 ) 

and evaluate the expectation values: 



f ^m+m^A'v) 2 • (62) 

k'X' 

Requiring again £(k) to be finite for k — > and integrable (see section 2) power counting reveals 
that the number of particles per mode possess a leading term N^ x ~ u>^ at low energies leading 
to a non-thermal spectrum e{uS) = 0{oj a ) also in this case. 
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With the aid of the completeness relation for the polarisation vectors e^ x and e^ = k/k : 

% ® % + E k\ ® Hx = 1 ( 63 ) 

A 

we arrive at 



f^m+u)?(i-{^) 2 ) • (64) 



For reasons of simplicity we consider radialsymmetric functions £ (r , t) — £ (r, £) and expand the 
Fourier transform £ into a Taylor-series: 

dt e iW * / d 3 r e ifef f (r, t) = / dt e icJ * £ A^„(i) (P) n = ]T M n {u) {k 2 ) n (65) 

^ ^ n=0 n=0 

assuming that the M n (t) which are related to the Mellin transform of £(r, t) exist, i.e. that 
£(r, t) is sufficiently localized. 

The zero term is proportional to the volume V of the disturbance. For example for a bubble 
described by e(r,t) =e oc + (l- e 00 )6( J R(t) - r) it follows: M (t) = ~(1 - e£)V(t) . 
Inserting Eq. (|65|) into Eq. (|6~4j) leads to: 



/•,.(. CO oo 

= f (1 - (^%a) 2 ) E E M* n {u % + uQ M m {u % + uQ (k + k>r n+m) ■ (66) 



n=0 m=0 



For the radialsymmetric case Nz x does not depend on A and eV , therefore we turn to Nk 



A 



N k = t -jrf (l + (%%0 2 E E A<>jE + -H»("S + u/ s ) (k + k'f {n+m) ■ (67) 



n=0 m=0 

fc' 



Evaluating the total radiated energy 



E= £ ujrNk -> / d 3 fc -LiV fc (6* 



(27T) 
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we arrive after some calculations at 

oo oo / / 7 \ 4+2n 



n=0 m=0 



4+2n \ / / 7 \ 4+2m 

a 




with a geometry-independent matrix C/ nm = Q nm {n + 171,600) (see appendix) 
The lowest (volume) term becomes for the bubble example 



If we assume that the dynamics of the disturbance can be described by a maximal characteristic 
length scale -R max and minimal time scale T min , this term is of order 0(R^ x T~7). Further terms 
contain additional factors of order 0(R^ aax T~£) and can be neglected for the case -R max /T miI1 <C 1. 



5 Velocity effect 

As one can see in section 2, the spectrum of the photons produced by the dynamical Casimir 
effect depends very strongly on the concrete form of the disturbance. Now we shall study some 
examples for the profile of the velocity j3(r,t): 

The first example corresponds to the radialsymmetric flow of an incompressible fluid around a 
bubble with the time dependent radius R(t): 

j3(r,t) = R(t) e r R 2 (t)/r 2 . (71) 

The Fourier transform f3 (k, to A remains finite in the limit k — > and is of order -R^ ax . Therefore 
the spectrum of the created photons behaves as uj a + 0(uj 5 ) for smooth and localized functions 
R(t) and the velocity contribution is for this case of the same order as the squeezing term. 
(The perturbation Hamiltonian for the squeezing effect is also of order -R^ ax due to the d 3 r- 
integration.) For the incompressible case the flow § fid A of the fluid is instantaneously constant 
even for |r | — > 00. If we assume a more realistic situation, where /3(r,t) yields relevant con- 
tributions only over a bounded volume of order -R^ ax , the velocity effect turns out to be of 
order and can therefore be neglected in those cases, where i? max is sufficiently small. (The 
^(i^) -behaviour results from the (i 3 r-integration that implies C*(-R^ ax ) and the assumption, 

— * 

that P(r,t) is still of order i? max .) 
For the second example: 

0(r, t) = R(t) e r (72) 



that is studied in Ref. [10], one obtains a Fourier transform that is divergent for k — * 0. This 
demonstrates the highly nonlocalized character of this profile which results in an increasing flow 
for |r I — > 00. Of course, for that case one could not provide the prediction e(uS) ~ oj a + 0(u> 5 ). 
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For a third example we study the situation of a rigidly moving medium with 

e(f, t) = Boo and 0(r, t) = (3(t) . (73) 

For that case Eq. (|34D for $ is easyly solvable: $ = (e^ 1 - + O0 2 ). Calculating the 
perturbation Hamiltonian H\ we find no photon production to first order in (3 (nonrelativistic 
limit): H\ |0) = 0{f3 2 ) . This implies that for the case of a purely time-varying velocity the 
effect will be of order N Rx = O0 4 ). 



6 Discussion 



6.1 Finite temperatures 

The response theory approach to quantum radiation allows for the incorporation of finite tem- 
perature effects. Instead of using the time evolution for the pure quantum state Eq. (|) one 
has to consider the equation of motion for the statistical operator p: 

and Eq. fll2[) becomes 



(^£4) = Tr I A) 7 " ex P (i J d 4 xHi(x] 



N% A r 



exp 



d x' Hi(x] 



(75) 



The zero-order term contains the undisturbed spectrum arising from p . The contribution of 
first order in H\ vanishes, if we assume a po related to the thermodymamic equilibrium: 



Po = p(t -> -00) = exp(-PH )/Z 



(76) 



The terms of second order contain the production of particles and the interaction between the 
thermal radiation field and the disturbance. 



6.2 Permeable media 

The approach presented above is generalizable to the case of media described by e(f, t), /i(f, t) 
and P(f,t), if one adds a term proportional to (// — 1) F^ V F V p u p with the dual tensor 
i<V= e^upaF^ to the Lagrangian. 

For pure magnetic disturbances (e = 1) the calculations proceed in the same way, if one replaces 
the potentials A and $ by the dual potentials A and T (see e.g. |T7|), defined via: 

l3 = VxA , H = A + Vr . (77) 

A possible application could be the case of magnetic phase transitions, e.g. in connection to 
superconductivity. 
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6.3 Contribution to sonoluminescence 



There have been various attempts to interprete the phenomenon of sonoluminescence p2f as 



an effect of quantum radiation fl23| , |9], 10 1. In order to estimate the contribution of QR it is 
important to determine the characteristic parameters (time- and lenght scales) involved in the 



dynamics of the bubble [£4|, Pqj . Calculations of the stationary Casimir energy (see Refs. p6 



and p7| |) cannot account for the dynamical character of the effect. Without an upper limit for 
the volume -R^ ax and the time interval T max of the disturbance and a frequency cutoff K c the 
radiated energy could diverge due to the space-time integration and the mode summation, if the 
perturbation functions are not localized and smooth. But if we assume a localized disturbance 
and a frequency cutoff due to the fact, that e differs only from 1 for wave numbers k < K c , 
an estimation of the upper limit for the order of magnitude is possible. We envisage Eq. (|62|) 
as an example for the general case of the quadratic response, assuming that the higher order 
terms are much smaller: 



N ix = $^m + M)mA>x'f= f |/A/dVH(M,^A',,,,') (78) 

k'X' k'X' 

with a function H that is of order 1. The order of magnitude of the space-time integrations and 
the mode sum can be estimated using i? max , T max , K c and the quantization volume V: 



J d A x -> 0(T max RlJ and ^ -> 0(VK 3 c ) . (79) 

k 

This leads to = ^(^L-RL^c/^) which depends still on the quantization volume V. 

But in the total energy as a physical observable V of course does not enter: 



E™*= f^ k Nf^ - 0{Rl a Tl,Kl) ■ (80) 
%\ 

Another result of this approach is the fact that the Casimir contribution to sonoluninescence 
should display a c<j 4 -spectrum for small u if we assume that the suppositions of section 2 are 
fulfilled. This provides a possibility to distinguish this effect from others (e.g. photon emission 
due to heating the gas in the bubble) which behave like u 3 . 



16 



Acknowledgement 

G. S. acknowledges support by BMBF, DFG and GSI. 

Appendix 

We have to introduce a complete set of real, orthonormal and transverse eigenfunctions / Q (r ) 
of the Laplace operator: 

V 2 A(r ) = -n 2 J a (r ) , (81) 

Va Vja(f) = , (82) 

Jd 3 rf a (r)f b (r) = 5(a,b) , (83) 

f f l a (r)m?>) = 5f(r-r>) . (84) 

a 

Accordingly, the field modes also satisfy the relation: 

J d 3 r (V x f a (f)) (V x f b (f)) = Q 2 a 5 ab . (85) 
Let us show how to calculate the canonical momenta p as well as their field representation 

— * — * 

II = Pafa from the Lagrangian: 

dL r jS dC r ,o dC dE 
dq a J dq a J dE dq a 

= J d 3 r D(f) (f a (f) + V X „(r, tj) = J d 3 r D(f)f a (f) . (86) 
For the same reason $ yields no direct contribution to , i.e., to V x H. 

W)=p a fa{r) = J d 3 r' f a (r) (f a (r')D(f')) 

= J d 3 r' e t 5f(f- r')Dj(f') = D ± (r) = D(f) (87) 

The transverse part F±(f,t) of an arbitrary vector field F(f,t) = F\\(r,t) + F±(f,t) is defined 



as 



V 2 F ± (f ) = -V x (V x F{ 
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As a complete and orthonormal set of real and transversal modes we use 

WO = /|cos(£f)% A (89) 



and 



fkX-(r) = \I^Mkmx (90) 



with polarization vectors eV, = e_^ A . The q are related to the creation/annihilation operators 
used in Eq. (^), e.g. for q^ x+ it follows: 

<% A+ = ^-==(af A + %a + + «-feA) • ( 91 ) 
To be complete we note here the explicit expression for the Matrix Q nm used in section 4: 



El "t" ) 2'((-l)'+ l) fv-J-4- 1 



(27r) 3 (8 + 2n + 2m)! ^ \ i ) v + l £ + 3 



n+m—£ 



Y, ( I (4 + £ + 2s)!(3-£-2s + 2n + 2m)! . (92) 



s=0 
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